Variable-range hopping in 2D quasi-lD electronic systems 
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A semi-phenomenological theory of variable-range hopping (VRH) is developed for two- 
dimensional (2D) quasi-one-dimensional (quasi-lD) systems such as arrays of quantum wires in 
the Wigner crystal regime. The theory follows the phenomenology of Efros, Mott and Shklovskii 
allied with microscopic arguments. We first derive the Coulomb gap in the single-particle density 
of states, g(s), where e is the energy of the charge excitation. We then derive the main exponential 
dependence of the electron conductivity in the linear (L), i.e. <r(T) ~ exp[— (Tl /T) lL ], and current 
in the non-linear (NL), i.e. j(£) ~ exp[— (£nl/£) 1nl ], response regimes (£ is the applied electric 
field). Due to the strong anisotropy of the system and its peculiar dielectric properties we show 
that unusual, with respect to known results, Coulomb gaps open followed by unusual VRH laws, 
i.e. with respect to the disorder-dependence of Tl and £nl and the values of 7l and 7jvl- 
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I. INTRODUCTION 

Issues related to transport in low-dimensional elec- 
tronic systems are challenging, the main difficulties be- 
ing to take into account of disorder or interactions and 
in some cases of the interplay between them. As known 
from Mott and Towse [lj, for strictly one-dimensional 
(ID) systems, any disorder leads to localized states. Such 
a statement implies that localization and transport prop- 
erties of ID Anderson insulators may be tackled per- 
turbatively by considering the limit of a weak disorder 
(so called Gaussian disorder where the density of im- 
purities N — > oo, their strength W — ► while NW 2 
is constant). For such a weak disorder, and in the 
absence of electron-electron and electron-phonon inter- 
actions, Berezinskii has confirmed the statement of 
Mott and showed that the a.c. conductivity is given by: 



' log lj. Subsequently, his approach has been 



extended by Gogolin, Mel'nikov and Rashba (GMR) 3] 
to the case where an electron-phonon coupling is present. 
They have shown that three-dimensional (3D) phonons 
provide a delocalization mechanism for the electrons at 
temperatures low enough that the scattering is mainly 
elastic, e.g. Tj„(T) ^> r e ; where Tj n (T) is the inelastic 
phonon scattering time and r e ; is the elastic scattering 
time on the static defects. This delocalization mecha- 
nism leads to a power-law hopping for the d.c. conduc- 
tivity, e.g. o~q(T) ~ T 3 , where the power comes from the 
phonon-scattering time. It should be noted that their 
arguments are not valid at the lowest temperatures; in 
particular, they hold only above To oc l/r e j, here and 
below: K = 1 unless specified. Subsequent studies aimed 
at exploring the effect of electron-electron interactions 
in ID Anderson insulators. In this respect, it has been 
shown Q that, in a Luttinger liquid, the Gaussian dis- 
order is strongly renormalized by interactions. For re- 
pulsive interactions, each impurity becomes effectively 
strong. The ID disordered interacting system is then 
equivalent [5j to an ensemble of weak links where impu- 
rities act as wire breakers. As a consequence, the power- 
law hopping laws acquire a non-universal exponent, e.g. 
interaction-parameter dependent. More recently 0, the 



low-temperature situation where the coupling to phonons 
is absent and electron-electron interactions dominate has 
been addressed. The modern notion of dephasing, due to 
these electron-electron interactions, has been considered 
as the delocalization mechanism of the electrons. It leads 
to a power-law hopping regime, reminiscent of the results 
of GMR, followed by a drastic suppression of the conduc- 
tivity Q- Less explored, from microscopic techniques, 
is the low-temperature regime: T <C To cx l/r e i, with 
both electron-electron and electron-phonon couplings. In 
this case, it is generally believed that the transport is of 
the variable-range hopping (VRH) type. In systems with 
3D phonons, that we shall be concerned with in the rest 
of this manuscript, semi-phenomenological arguments by 
Mott suggest that the VRH laws are characterized by 
stretched exponentials and read: 



a(THao(T)exp[-(T /m 



(1) 



where 7 < 1. In ID, To oc where the single-particle 

density of states v oc 1/vf and the localization length 
£ cx vpT e i. Hence, To oc l/r e /, and the exponential de- 
pendence manifests at T <C To oc l/r e j. In the oppo- 
site case, the temperature dependence of the conductivity 
arises mainly from the pre- factor: 00 (T), which depends 
on the electron-phonon coupling, and leads to power-law 
hopping, see Fig. ^for a summary of regimes. The suc- 
cess of Mott's arguments came from their wide experi- 
mental confirmation for a great variety of disordered sys- 
tems, i.e. it has been found experimentally that 7 « 0.25 
in isotropic 3D and 0.33 in isotropic 2D systems, in ac- 
cordance with the theory. Later, Efros and Shklovskii 
(ES), see Ref. Q for a review, extended these arguments 
to the Mott-Anderson insulators where, besides strong 
disorder, the long-range Coulomb interaction is present. 
Once again, their results were confirmed by a great vari- 
ety of experiments in the field of doped semiconductors 
where the exponent 7 ss 0.5 in all dimensions^ as pre- 
dicted by ES. 

In the present study, we focus on the low-temperature 
transport properties of strongly disordered and inter- 
acting two-dimensional (2D) quasi-lD electron systems. 
The latter consist of a 2D periodic lattice of parallel 
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FIG. 1: d.c. conductivity as a function of temperature for 
a quasi-one dimensional system coupled to three-dimensional 
phonons. In the range of temperatures displayed, phonons 
provide a derealization mechanism for the charge excitations 
of the system. Power-law hopping (PLH) dominates at tem- 
peratures T 3> To and variable-range hopping (VRH) domi- 
nates at T <C To (Tb c< 1/T e i where r e ; is the elastic scattering 
time on the static defects). 

wires and therefore display markedly anisotropic prop- 
erties with relations to the ID world, see Refs. [TlJ and 
[l2| for reviews on ID physics. We assume that im- 
purities are point-like and act as wire breakers. Inter- 
actions are long-ranged and lead to a large parameter: 
fs = Uh/£f, corresponding to the ratio of the Coulomb 
energy-scale: Uh = e 2 / na (k is the dielectric constant 
of the host lattice in which the 2D system is embedded), 
to the kinetic energy scale given by the Fermi energy: 
Ef- In the field of charge-density waves (CDW)i£, such 
quasi-lD Mott insulators are referred to as Akp CDWs. 
This Akp modulation corresponds to a space periodicity 
of the charge-density along the wires, X^kp = 2n /Akp 1 
equal to the average distance between electrons a; such 
systems therefore display a strongly correlated state of 
the Wigner crystal type. For our purposes, such systems 
include man-made atomic or molecular wire arrays em- 
bedded in a s emi-co nducting matrix [l5j and mono-layers 
of CDWs [l6lll7ll8j . For these quasi-lD Mott- Anderson 
insulators, where the interplay between disorder and in- 
teractions is quite non-trivial and out of the reach of any 
perturbative scheme, we will naturalljiiS follow the semi- 
phenomenological route of Mott and Efros and Shklovskii 
in deriving the low-T (i.e. T <C l/r e j)Si transport prop- 
erties. Our motivation resides in the existing experimen- 
tal literature on quasi- ID systems where VRH la ws hav e 
been reported quite extensively, e.g. in polymers |22l23| |. 

Eq. corresponds to linear variable-range hopping 
laws such that the current j is proportional to the applied 
electric field £, j = a(T)£. Upon increasing the applied 
electric field, such laws cross-over to non-linear VRH laws 
which read: 

j(£)~exp[-(£ /£) 7 L (2) 

where the exponent 7 < 1. In the frame of disor- 
dered semiconductors in high electric fields such laws, 
with an exponent 7 = 1/4, are known from Mott and 



Shklovskii. In quasi-lD systems they are observed, e.g. 
for bronzes [24[. These laws hold in the low-T limit, 
T <C 1/Vei, up to the threshold field for global slid- 
ing of the pinned electronic structure, see Refs. [25j and 
[26j for reviews on pinning. Of course, for the exponen- 
tial character to manifest the electric field should not 
be too high, e.g. in ID: £ <C £q oc 1/evFT 2 ^ For 
a given linear VRH law, it will however manifest for: 
£ > £ c where the crossover field reads, in any dimension: 
£ c = T£q/Tq oc T/e£, where e is the unit charge and £ 
is the localization length. In summary, we shall focus on 
the following two regimes: the linear response regime of 
Eq. JTJ: 

T < T « 1/rei and £ < T/e£, (3) 

and the non-linear response regime of Eq. : 

T < T « 1/rei and T/e£ « S < So- (4) 

Some of the recent theoretical literature on VRH has 
been devoted to ID and 3D quasi-lD systems in the lin- 
ear regime, cf. Rcf. 14] . To the knowledge of the author 
a theory of electron transport in 2D quasi- ID systems 
has not been developed in either the linear or non-linear 
regimes. Our specific task will be to compute the main 
exponential dependence of the electron current as a func- 
tion of temperature, electric field and impurity concen- 
tration. We will closely follow the arguments of Ref. n4 . 

Our results display two characteristic features: a non- 
monotonous dependence of the current as a function of 
disorder and a highly non-universal exponent 7, i.e. im- 
plicitly interaction- and disorder-dependent. Both state- 
ments will be proved in the following sections. They are 
generic of incommensurate quasi-lD Mott-Anderson in- 
sulators and are absent in the usual doped semiconduc- 
tors or even in the other (hypothetical) collective struc- 
ture of interest: the pinned Wigner crystal. The reader 
interested more in our results than the details of the 
derivations may refer directly to Sec. IVIII 

The theory will be developed as follows. In Sec. [H] 
we describe the 4/cf quasi-lD systems in a semi-classical 
way, which is justified due to the large parameter r s . 
This will unable us to determine the dielectric properties 
of the system. With the help of these results and the ES 
phenomenology we derive, in Sec. 11111 the single-particle 
density of states (DOS) of charge excitations and show 
how it is affected by the long-range interaction potential 
among charge excitations (Coulomb gap). In Sec. II VI we 
focus on determining the localization length with the help 
of microscopic models related to single-impurity tunnel- 
ing of the charge excitations. In Sees . M and I VII the VRH 
laws, in the linear and non-linear response regimes, re- 
spectively, are derived with the help of the Efros, Mott 
and Shklovskii phenomenology and the results of Sec. lIIII 
and Sec. II VI In Sec. IVIll the conclusion is given. 
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II. DIELECTRIC PROPERTIES 

In this Section we determine the dielectric properties 
of disordered 2D quasi-lD systems. The latter will be 
used in subsequent sections to determine the VRH laws. 

1. The model 

Following the Introduction, we introduce our basic 
model of a strongly pinned 2D quasi-lD system in a 
4k f — CDW state. Along each wire, such a quasi-lD sys- 
tem is characterized by a modulation of the density: 

p(x) = po + pi cos(Qx + ip) + (l/2%)d x (p, 

where po is the unperturbed density, Q = 4fcp is the 
modulation of the wave, ip is the phase of the (assumed 
incommensurate) CDW and the last term describes long- 
distance deformations. In electronic systems, a 4kp mod- 
ulation implies that the wavelength of the CDW is of the 
order of the average spacing a between electrons. It is 
realized in systems where the long-range Coulomb inter- 
action is present, see Ref. [llj for a review on such results 
in the ID case. In quasi-lD systems, where the Coulomb 
interaction may be screened by neighboring chains, the 
4fcp-CDW is also detected [l7|. This reflects the large r s 
nature of such systems. It implies that one may neglect 
quantum fluctuations and describe these CDWs semi- 
classically with the help of the phase-field ip. This CDW 
phase, ip, is related to the order-parameter describing the 
condensate which reads: A(r) = |A| exp[iip(r)], where r 
is a two-dimensional coordinate, and |A| the amplitude 
which is assumed to be frozen as we consider low temper- 
atures. Moreover, the system has charge invariance, i.e. 
ip — > ip + 2tt, implying the existence of a Frohlich mode, 
i.e. the sliding of the electron crystal, in accordance with 
the fact that the system is incommensurate. 

We focus first, in the frame of this semi-classical the- 
ory, on static properties of the system: the determination 
of the structure of the electron system in the presence of 
a single strong impurity and of it's charge excitations as 
well as the electrostatic potential between such excita- 
tions [the dynamics will be considered later, i.e. the ef- 
fect of quantum fluctuations, when dealing with the tun- 
neling of these charge carriers through the impurities]. 
Adding then a single strong impurity at the origin of the 
2D system, the Hamiltonian consists of three parts: 

H = H + H D +H I , (5) 

where Hq is the elastic part and reads: 

Ho = J d 2 r | [(d x ip) 2 + a(d y ip)\ (6) 

where Y — Y x is the bulk-modulus along the chains and 
the anisotropy parameter reads: a — Y y /Y x <C 1 (x— is 
the direction along the wires and y— is the direction per- 
pendicular to them) . In Eq. JfJJl , the phase field ip may be 



interpreted as a scalar displacement field along the wires, 
u x , as in conventional elasticity theory. The displacement 
and phase fields are then related by: ip = —(2ir/a)u x . It 
follows naturally that the first term, in Eq. © , describes 
the compression energy along each wire, Y being the bulk 
modulus and the second term corresponds to shear elas- 
ticity due to inter-wire interactions, Y y being the shear 
modulus. It should be noticed that the inter-wire inter- 
action term is defined at the level of the order parameter 
with the help of: A*A m + c. c. = |A| 2 cos(y„ — ip m ), 
where n and m index neighboring wires. Shear elasticity 
is then derived by expanding the cosine, at low temper- 
atures, and going to the continuum limit perpendicular 
to the wires. 

The second term in Eq. (J5J describes the effect of a 
strong impurity at the origin and reads: 

Hd = J d 2 r [W f d x pS(r) - W b eos{ip)5(r)], (7) 

where the first term corresponds to the forward scatter- 
ing on the impurity (i.e. the coupling of the point- like 
impurity potential to the long-distance part, cx d x ip, of 
the CDW density) and has a strength Wf. The second 
term corresponds to the backscattering on the impurity 
(i.e. the coupling of the point-like impurity potential to 
the oscillating, 4kp, part of the CDW density) with a 
strength Wb- 

Finally, the third term in Eq. J5J contains the long- 
range Coulomb field and reads: 

Hi = J d 2 r [U5(r) + ±d x ipU(r)\ d 3 r(VC/) 2 ],(8) 

where U is the Coulomb field. The first term corresponds 
to a point-like test charge for the Coulomb field U. The 
second to the coupling of the long-range Coulomb poten- 
tial to the long-distance part of the density (<x d x ip) and 
the last term corresponds to the energy of the Coulomb 
field. It should also be noticed that, in the model of 
Eq. (JSJ), the Coulomb interaction is the real one, i.e. the 
three-dimensional one. 

There is a basic non-trivial length-scale in Eq. © 
which is proper to quasi-lD systems and will bring signif- 
icant opportunities to go beyond strictly ID physics all 
along this manuscript. This is the length of a 2tt— soliton: 
l s , and emerges due to the fact that in quasi- ID systems 
shear and elasticity are coupled 27 . This can be under- 
stood with the help of the following arguments. An impu- 
rity will enforce a deformation of the CDW in its vicinity. 
In the absence of long-range Coulomb interaction, a de- 
formation Sip along a distance Sx along the defected wire, 
has an elastic energy: E — Sx b 2 [YSip 2 /Sx 2 + Yj_5ip 2 /b 2 ], 
where we have assumed that the distance between neigh- 
boring wires is b. Minimizing this elastic energy with 
respect to Sx and defining the optimal Sx as the length 
l s , we find that: 

Is = -JU ( 9 ) 
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In the case where Sip = 2n, the deformation becomes the 
plastic one, as a whole period of the density-wave is af- 
fected by the impurity. Hence the name of soliton length, 
see Refs. [2^ and [2^ for reviews on soliton physics in 
condensed matter. Notice also that for decoupled wires, 
i.e. strictly ID systems where a — > 0, l s is infinite. This 
implies that the whole wire adjusts to the deformation 
at the origin. In the case of < a < 1, inter- wire inter- 
actions enforce the same phase [2tt] between neighboring 
wires, beyond the length on both sides of the impurity. 
In what follows, the length l s will be used to distinguish 
between the low- impurity concentration limit, where the 
average distance between impurities along a wire, Z, is 
larger than l s , I ^> l s and the opposite large-impurity 
concentration limit, I <C l s (the density of impurities is 
given by: N = 1/lb 2 ). 



2. The case of small impurity concentration 



in the quasi-lD system. As can be seen from Eqs. Ill all 
and 1|12|) . in 2D, the effect of the Coulomb interaction 
is to shift q 2 to g 2 £2(q) = q 2 + q 2 /rrjq. This shift can 
be included in the bulk modulus and corresponds to the 
incompressibility of the crystal as a whole. As we are 
interested in physical properties at large distances along 
the chains q x <C \faq y or 1 > y/yfa, we see that the 
Coulomb interaction brings into play the following sector: 



q 2 x <^ar D ql y < {ar D x 2 ) 1/3 , (13) 

which plays a dominant role as we show now. 
The inverse Fourier transform of Eq. i|llfl yields: 



ip(r) = 
f/(r) = 



—b sgn(x) 

6ir(y?arD\x\) 2 / 3 
2 

1 — exp 



e 

k\x 



3r ( f>- 2 r ( f,Vp™ 

y/ar D \x\ 
Tfi/ 2 



(14a) 
(14b) 



We focus on the low-impurity concentration case, I — 
1 /Nb 2 ^> l s (recall that I is the average distance between 
impurities along a wire), cf. the one-impurity model of 
Eq. ©: 



1 

T7b 



b 2 



(10) 



Each impurity is strong and the backscattering term in 
the Hamiltonian describes the pinning of the system at 
distances of the order of the average distance between 
impurities, I. On the other hand the forward scattering 
term describes softer deformations of the system away 
from the impurities [cf. the dipole distortions that we 
will introduce below]. We focus on the lattes^ and set: 
Wb = 0, assuming that the system is strongly pinned. By 
varying the functional of Eq. JSJ with respect to the phase 
field and the Coulomb field, 5H/5<p = and SH/SU = 
and going to Fourier space, one gets a system of two 
Poisson equations where both fields screen each other: 



</>(q) = ~ 



iqx 



<f £2(q) 



2ne 2 
q?2 (q) 



Wj_ 

Y 



b 

ruq 



(11a) 



(lib) 



where q 2 = q 2 + aq 2 , q 2 = ql + q% and r^ 1 = 2ne 2 /Yb 2 w 
6 _1 is the inverse screening length in the metallic phase. 
The 2D dielectric constant with the help of which the 
long-range Coulomb and elastic fields screen each-other 
reads: 



£ 2(q) 



1 



r D qq 2 



(12) 



which displays the anisotropic metallic-like (in the sense 
of sliding of the collective electron structure) screening 



where T(x) and T(x,y) are the complete and incom- 
plete gamma functions, respectively. The phase field 
of Eq. (|14a(l is an odd function of the coordinate along 
the wire and at long distances along the wire, i.e. for 
y <C (arrjx 2 ) 1 / 3 , reads: ip(r) oc -(l s /x) 2 / 3 , cf. Eq. I|15a|) 
below. This implies that, at such distances, the impurity 
is surrounded by dipole distortions of length Z s . Cru- 
cially, we identify these dipole distortions as the charge 
excitations of the system. 

These charge excitations interact with a Coulomb po- 
tential which deviates from the usual 3D one, due to 
the anisotropic screening, cf. Eq. jug. More gen- 
erally, Eq. Ijllbj) shows that the electrostatic potential 
changes sign along the cone: y — ±y/a\x\/ (Wf /Yb — 1). 
This is related to the dipole nature of the phase defor- 
mations. In particular, outside this cone, the potential 
reads: U cx — Wfrrj/r 3 , so that charge deformations of 
the same sign attract each-other via a dipole potential. 
On the other hand, within this cone (which is the sec- 
tor of validity of Eq. 114b|l ) the potential is repulsive 
and independent of Wf. This implies that, if they were 
not bound to the impurities (which they actually origi- 
nate from), such charge excitations would form domain 
walls (our arguments have their roots in Ref. |32j). In 
the following we return on the sector within the cone: 
y < ±^/a\x\/(Wf/Yb— 1), where the potential is repul- 
sive and Eqs. Ijl4|l holds. 

We summarize this sub-section by giving the asymp- 
totic expressions of Eqs. (|14f> . on one hand close to the 
chains, i.e. in the sector of Eq. ljT3|) . y (arrjx 2 ) 1 ^ 3 , 
where it reads : 



ip(r) = -sgn(x) 



2irx 



2/3 



U(r) 



k\x\ ' 



(15a) 
(15b) 



and on the other hand further away from the chains, i.e. 
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y 3> (arox 2 ) 1 / 3 , where it reads: 



<p(r) = - 



b sgn(x) 
27r|y| : 



U(r) 



3/2 



(16a) 
(16b) 



In the limit of vanishing inter-wire coupling: a — > 0, the 
sector defined by Eq. I|13(l vanishes. This corresponds 
effectively to a crossover from a low impurity density 
regime, N <C 7V S = l/l s b 2 — y/a/b 3 , to a large impu- 
rity density regime, iV 3> iV s , that we consider next. 



3. The case of large impurity concentration 

Formally, this CctSCj clS defined in the previous section, 
requires that I <C l s (recall that I is the average distance 
between impurities along a wire) or, in terms of impurity 
concentration, that: 



iV> N s 



IP' 



(17) 



where a gives the dimensionless strength of inter-chain 
interactions. Large impurity concentration is therefore 
equivalent to vanishing inter-chain couplings. The sys- 
tem is then equivalent to an ensemble of metallic (in 
the sense of sliding) segments along the wires, of aver- 
age length I. These segments are decoupled elastically 
but still coupled by the long-range Coulomb potential. 
In the literature, such a regime is sometimes referred to 
as a model of interrupted metallic strands^!. Following 
the previous paragraph, our goal here is to determine the 
dielectric properties of such a phase. 

From Eqs. Ijllfl . the electrostatic potential of a charge 
carrier at the origin of a pure system (Wj — Wb = 0) 
reads: 



2ne 2 
?£2(q)' 



(18) 



where the dielectric constant is given by Eq. H12|) and 
Q 2 = Qx + Qy- ln the previous paragraph we have fo- 
cused on the effect of deformations due to the forward 
scattering term (Wf) in a system of dilute strong pin- 
ning centers. Increasing the impurity concentration we 
need to take into account the average effect of these pin- 
ning center. This is done according to the prescription: 
q 2 — > q 2 + L x 2 in Eq. I|12f> for the dielectric constant (re- 
call that q 2 = q 2 + aq 2 arises from the elastic part of the 
energy) where L x oc I is the pinning len gth due to the 
backscattering on the impurity, see Refs. l25l and [26 for 
reviews on pinning. This amounts to introduce the effect 
of backscattering on the impurity as a commensurability 
term explicitly breaking the translational invariance of 
the system. The dielectric constant therefore becomes: 



£ 2(q) 



1.7- 



roq(q 2 + L x ) 



(19) 



In the clean limit, L x oc I ^> l s , we may assume that 
L x /l s — * oo, and recover the results of the previous sub- 
section. In the dirty limit, L x oc I <c l s , we may assume 
that l s /L x — > oo so that a — + and q 2 = q 2 + aq 2 — > q 2 . 
In this regime, the Coulomb potential reads: 



U(r) 



2Tre 2 



(20) 



The asymptotics of Eq. 1121 III then read: 



U (r) = —r 1-exp - 




U(r) 
U(r) 



^JybK x /K / 
L x < x < L 2 x /2b, y<^L 2 x /2b, (21a) 

-, L x <Ci< L x /2b, y > L 2 x /2b, (21b) 



ny 

e 2 

— , r»L2/26, 
Kr 



(21c) 



where r = y/x 2 +y 2 and k x is the longitudinal (along 
the wires) dielectric constant: 



k x = n(L x /r D y 



(22) 



up to a numerical constant of the order of unity, with: 
ro ~ &, where b is the inter- wire distance and the pinning 
length along the wires: L x w /, is of the order of the 
average impurity-distance, I, along each wire. Notice that 
the longitudinal dielectric constant, Eq. (|22|l . is inversely 
proportional to the square of the impurity concentration 
in the system, N = 1/lb. This longitudinal dielectric 
constant is large, i.e. larger than the dielectric constant 
k of the host in which the 2D system is embedded. This is 
important and explains the richness of Eqs. J2J. Indeed, 
Eq. I|21a|l shows that close to the chains (y <C I 2 /2b) and 
for I <C x <C I 2 /2b, field lines prefer to remain in the 2D 
layer where they are screened by the large longitudinal 
dielectric constant. This sector is peculiar to the quasi- 
1D system. At larger distances, x I 2 /2b, the lines 
escape to the external media where they are screened 
only by k and the potential is fully long-ranged. This 
result is derived by other means in the Appendix. 



III. SINGLE-PARTICLE DENSITY-OF-STATES 
AND THE COULOMB GAP 

We now turn to the determination of the single-particle 
density of states, 3(e), of the localized charge excitation 
of energy e. This density of states is defined as the prob- 
ability density for the excitation to have energy e. In the 
case where e = £j, the corresponding wave function is 
given by: 



Mr) 



exp[-— - — ] 



(23) 



where £ x is the longitudinal, i.e. along the wires, local- 
ization length and t; y is the transverse one. Moreover, 
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Xi and yi locate the position of the localized charge ex- 
citation of energy Ej. Notice that, in usual disordered 
semiconductors, V is the single-particle wave function of 
the electron (impurities are hydrogenoid atoms) and the 
localization length corresponds to the Bohr radius of this 
trivial charge excitation. In quasi- ID systems ip is the 
single-particle wave function of deformations of the elec- 
tronic system, the distortions considered in the previous 
section. 

An ensemble of such localized states with a density 
g(s), forms a disorder or impurity band. Moreover, keep- 
ing in mind that the system is a Mott insulator, the local- 
ized states of energy e belong the lowest Hubbard band 
of the system. It is on this lowest "Hubbard impurity- 
band" that we shall focus now, keeping in mind that the 
upper one, distant by the Hubbard: Uh ~ e 2 /na, is un- 
reachable at the low energies we consider: T -^Uh- 

We further assume that this lowest Hubbard impurity 
band is partially filled. As a result, the low-T transport 
we consider is the hopping of charge excitations within 
this lowest Hubbard impurity band. One may formally 
define a Fermi energy for this band of localized states 
(ep = 0) and a non-zero DOS of charge excitation at 
this Fermi level: v = g(0) ^ 0. As far as transport is 
concerned, we shall show below that the non-zero i.e. 
the existence of gapless charge excitations, implies that 
there is a low-temperature hopping conductivity of the 
variable-range type. 

In the absence of the long-range Coulomb g(e) is con- 
stant, i.e. g(e) ~ is, within the disorder band width. 
We now turn on the long-range Coulomb interaction and 
follow Efros and Shklovskii to determine the influence of 
this interaction on g(e). We provide here some details 
for the reader which is not familiar with the ES argu- 
ments. While hopping from the localized state i of en- 
ergy Si below the effective Fermi energy of the impurity- 
band to the state j of energy ej above, the energy of the 
charge excitation will vary by: SE = Ej — e 2 /nr^ — E{. In 
the latter expression, the long-range Coulomb interaction 
— e 2 / ' KTij between the charge excitation at j and the hole 
it has left at i has been taken into account. By construc- 
tion: SE > 0, which implies a depletion of states around 
the Fermi energy: r« > e 2 /n\£j — £j|, because of the 
long-range nature of the Coulomb interaction. In d di- 
mensions, such states have a spatial density of n = l/r d , 
where r = rij. This implies that: n(e) < (K,\e\/e 2 ) d , 
where e = Sj — £j. The corresponding single-particle 
density of states (g(e) = dn(e)/de) then reads: 

g(e) = C(^) d \e\ d -\ 



affects the transport as will be shown below. Moreover, 
the Coulomb gap of Eq. I|24|) depends on dimensional- 
ity^: quadratic in the energy of the charge excitation 
in 3D (oc e 2 ) and linear in 2D (ex and not on the 
impurity concentration. The situation is much richer in 
quasi- ID systems due to their non-trivial dielectric prop- 
erties. 



4- The case of large impurity concentration 

We apply the Efros- Shklovskii arguments to the case 
of the anisotropic Coulomb interaction characteristic of 
quasi- ID systems starting with the large-impurity con- 
centration regime: N 3> N s . The Coulomb interaction is 
given by Eqs. 1|21[1 . For each sector we fix the potential 
U = e, determine the equipotentials x(s) and y(e) and 
substitute them in the density, n(e) in order to derive the 
Coulomb gap: 



n(e) 



1 



x(e)y{e) 



9(e) 



dn(e) 
de 



We first focus on the large- distance sector of Eq. I|21c|l . 
At such distances, x 3> I 2 /2b, the potential is isotropic. 
Therefore the usual ES arguments apply, x(e) = y(e) = 
e 2 / ke and n(e) = K 2 e 2 /e 4 . The Coulomb gap is therefore 
linear in the energy of the charge excitation: 



gEs{e) = C -j |e|, £ < £i, 



(25) 



where Co is a numerical coefficient of the order of unity, 
£i is a crossover to another Coulomb-gap shape (as ex- 
plained shortly). As in the usual ES law, Eq. 12511 is 
independent on the impurity concentration (Co and n 
are constants). 

We then focus on the shorter-distance sector of 
Eq. (|21aJ) . At such distances, I <C x < I 2 /2b and 
y <C I 2 /2b (recall that L x cx /), the potential reads: 



U(r) 
U(r) 



2/< 



bn x 



(26a) 



Ky/ybK x /K 



y»^. (26b) 



which follows from Eq. (|21a|l . Close to the chains, 
y <C x 2 K/bn x , and fixing U = e yields the following 
equipotentials: x(e) — e 2 / ke and y(e) — x{e) 2 k/ an x . 
This yields a density of localized states: 



(24) 



n(e) 



bK 2 K x e 3 



where C is a numerical coefficient of the order of unity, 
depending on the dimensionality of the systemifl, e.g. see 
Ref. |3^ for the 2D case we are interested in. A crucial 
feature of Eq. I|24|l is that the DOS vanishes only at the 
Fermi energy. Hence, there is still no Coulomb block- 
ade within the impurity band because of the long-range 
Coulomb. There is however a soft Coulomb-gap which 



and therefore the Coulomb gap: 

/ \ /-» bn 2 K x £ 2 
gi(e) = C x — -„ — , 



£i<£<Ai, (27) 



where C\ is a numerical coefficient of the order of unity. 
This Coulomb gap is quadratic for a 2D system and de- 
pends on disorder through k. x , cf. Eq. (|22|) . Both features 



FIG. 2: Schematic plot of the DOS of charge excitations 
of the 2D anisotropic system with a large concentration of 
impurities, N ~3> N s . The linear Coulomb-gap at low ener- 
gies [Eq. 1251 1 is succeeded by the quadratic Coulomb-gap 
[Eq. 1271 1 and finally by the constant DOS of the disorder 
band, v. 



FIG. 3: Schematic plot of the DOS of the 2D anisotropic 
system with a small concentration of impurities, TV <C N s . 
The linear Coulomb-gap at low energies [Eq. 1251 1 is succeeded 
by the quadratic Coulomb-gap [Eq. 1271 1. then by the "2/3- 
" Coulomb-gap [Eq. 1311 1 and finally by the constant DOS of 
the disorder band, v. 



are unusual with respect to known results, cf. Eq. (|24|l . 
In Eq. 1)2 7[). Ai is the width of the Coulomb gap and e-y 
a crossover energy from Eq. I|25|l to Eq. (|27|l . The upper- 
bound Ai to the energy-dependence of Eq. I|27|) origi- 
nates from the fact that this Coulomb gap is due to the 
short-distance part of the Coulomb potential Eq. (|21a|l . 
At shorter-distances, hence higher energies, Eq. i|27|) 
crosses-over to the constant DOS, v, of the disorder band. 
Equating <?i(Ai) to v yields: 



Ai = D x [e & v/bn 2 K x ] 



1/2 



(28) 



where D\ is a numerical coefficient of the order of unity. 
Eq. (|28|) corresponds to the Coulomb-gap width in the 
large impurity regime. On the other hand, the crossover 
energy E\ is obtained by equating gEsi^i) to gi(ei) 
which, from Eqs. and J23), yields: 



ei=JVi 



e 

6k, 



(29) 



where 7Vi is a numerical coefficient of the order of unity. 
Notice that both Eqs. (|28|l and (|29|l have an impurity- 
dependence through the longitudinal dielectric constant 
k x which is given by Eq. (|22|l . In the large impurity con- 
centration case, the total Coulomb gap shape is plotted 
on Fig.HJ 



5. The case of small impurity concentration 

We turn on to the small-impurity concentration case: 
/V <C N s . The ES Coulomb gaps found in the case of 
large impurity concentration, N N s , are still valid in 
the present case. Hence, at distances x larger than l 2 /2b, 
interactions open a linear Coulomb gap similar to the 
one of Eq. l(23|) . At smaller distances: ! « i « I 2 /2b, 
corresponding to larger energies, this linear Coulomb 
gap crosses over, at e\, to the quadratic one defined by 
Eq. I|27|l . Going to higher energies corresponds to dis- 
tances shorter than the pinning length, L x oc I. If impuri- 
ties are sufficiently diluted, i.e. I is large, we may face the 



situation where a charge excitation experiences the po- 
tential within the pinning area L x L y — \fol 2 , i.e. the po- 
tential of the equivalent pure system given by Eqs. (|15b|) 
and JT6EJ): 



k\x\ 



U(r) = - 

K 



ar D 



-1 1/2 



y < [x 2 ar D ] 1/3 (30a) 



y> [rrWc] 1 / 3 . (30b) 



Close to the chains, y <C [a^ar^] 1 / 3 , and fixing U = e 
yields the following equipotentials: x{e) = e 2 /ne and 
y(e) — [x 2 (e)arD] 1 ^ 3 ■ This yields a density of localized 
states: 



1 / KS \ 



(ar D y 

and therefore the Coulomb gap: 

£ 2/3 

52(£) = ° 2 (ar D )V3( e 2/ K )5/3' 



£ 2 <£<A 2 , (31) 



where Ci is a numerical coefficient of the order of 
unity. In Eq. I|31|) . the new crossover energy, e 2j de- 
fined as the crossover energy between the quadratic and 
2/3— Coulomb gaps \g\{ei) = 92(^2)} reads: 



e 2 


'1 /K x \ 3 ' 


Kb 


a V k J 



"I 1/4 



(32) 



where N2 is a numerical coefficient of the order of unity. 
Finally, A 2 corresponds to the Coulomb-gap width in 
the small- impurity concentration case [g 2 (A 2 ) = v\ and 
reads: 



A 2 = £ 2 [(e 2 / K )^ 3 ar D ] 



1 1/2 



(33) 



where £> 2 is a numerical coefficient of the order of unity. 
These results are summarized on Fig. [3] which displays 
the total DOS of the system in the small impurity con- 
centration regime. 
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IV. LOCALIZATION LENGTH AND 
TUNNELING 

We now focus on the localization length associated 
with the single-particle wave function of Eq. J23J). The 
transverse localization length is given by: £ y < b, where 
6 is the inter-chain distance. In the following we focus on 
£ x . As we have already said at the level of Eq. (J2HJ), in 
the case of doped semiconductors, the charge excitations 
are simply the electrons and the corresponding localiza- 
tion length is the Bohr radius of the hydrogenoid wave 
function. The later has no impurity dependence. On the 
other hand, the situation is much less trivial for quasi-lD 
systems, where excitations are non-linear deformations 
of the electronic system. These excitations are extended 
and one may think that they are localized over the length 
£ x = l s along the wires, in the clean limit (I 3> l s ) and 
over the distance £ x = I between impurities in the dirty 
limit (I <C l s )- This intuitive argument is correct but it 
turns out that the Fohlich mode as well as the long-range 
Coulomb interaction also influence, in a non-trivial way, 
this length, as we show now. The reader which does not 
wish to go through our microscopic arguments, at least 
on the first reading, may skip to the next section, where 
transport laws are derived and it is shown how the lo- 
calization lengths enter expressions for the conductivity 
and current. 



6. The model 

It is crucial to notice that the localization length is of 
purely quantum origin, contrary to the pinning length 
(Fukuyama-Lee length or Larkin length depending on 
the context) which is defined at the classical level. The 
derivation of the localization length requires to consider 
the tunneling of the charge excitations through impu- 
rities. The most convenient way of dealing which such 
processes, which are forbidden at the classical level, is 
to extend the semi-classical approach of previous sec- 
tions to the Euclidean space by using imaginary time 
dynamics. In doing so, we consider the tunneling of a 
charge excitation along the distance x of a given wire. 
The average distance between impurities along this wire 
is /. Hence, the charge has to tunnel through x/l impu- 
rities. Neglecting interferences between these impurities, 
their effect, as a first approximation, is additive. The 
total action for tunneling along the distance x is there- 
fore given by: S op t(x) = s op tx/l, where s op t is the action 
to tunnel through a single impurity. The total action 
S pt enters the quantum probability to reach the distant 
point x: |?/>(x)| 2 ~ exp(— S op t(x)), where ip is the one- 
(quasi-)particle wave function of Eq. (|23|l . Therefore: 
x/£ x = S op t = s ptx/l and the localization length reads: 



This result is valid for the large impurity concentra- 
tion case. When the concentration of impurity is small 
charge excitations acquire the length l s and the localiza- 
tion length becomes: 



opt, 



I > L. 



(35) 



In our approximation of neglecting interferences between 
various impurities, which is reasonable in the strong or 
individual pinning regime we consider, the localization 
length depends only on a one-impurity tunneling action. 
Subsequent calculations will therefore focus on a deriva- 
tion of this one-impurity optimal action s opt . 
The basic action that we shall consider reads: 



,h I »lr ! ^W) 2 + - [(d xV f + a {d^f 

(36) 



-W b cos(<p)S(r)}, 



where now ip = (p(r,r), C is the CDW stiffness, the 
last term corresponds to the backscattering on the im- 
purity and h — 1 unless specified explicitly At this point 
we came up to a systematic way of deriving the local- 
ization length with the help of Eqs. 1)340. I|35|l and the 
well-defined model of Eq. I|36l) . Unfortunately, it is not 
possible to determine analytically the exact non-trivial 
time-dependent solutions (instantons) of the above non- 
linear partial differential equation. The physical process 
of tunneling may however be understood as a two-stage 
process [34j. The first stage (small times) corresponds to 
the local tunneling of the CDW at the impurity position 
by one period, a. This stage is described by the backscat- 
tering term which is local in space. We now assume that 
the non-trivial dynamics of the phase at the impurity are 
essential only at large times (large and small times will 
be defined below). As a result the jump of the phase, 
at the impurity, during the first stage is described by the 
Ansatz: f(r) = 2tt9(t), where 9{t) is the Heaviside func- 
tion. The action describing the first stage is then linear 
in time and in the impurity strength: 



si = W b l3. 



(37) 



The second stage (large times), corresponds to the ad- 
justment of the crystal, at large distances from the im- 
purity, to the jump of the phase at the impurity. It is 
described by assuming that the local effect of the impu- 
rity is irrelevant, i.e. that one can set Wb — 0. Then, 
the corresponding effective action may be derived from 
Eq. H36|) (with W& = 0) by tracing out all remaining gap- 
less modes away from the impurity position. This leads 
to an effective action for the phase at the impurity: 



s 2 = ^pYl i^i 2j ^' 



(38) 



£x = 1 1 S opt , I < h 



where the spectral function J u contains the effect of the 
(34) infinite number of gapless modes which have been traced 
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out. Without the Coulomb interaction this spectral func- 
tion reads: 



dq 



1 



(2tt) 2 Clu 2 + Yq 2 ' 



(39) 



where q = q x + aq 



7. The case of large impurity concentration 

In this case wires are elastically weakly coupled, a — ► 
0. This case is worth examining first because it relates 
our approach to known results in the literature on ID 
disordered interacting systems. Eq. l(3l^l leads to J w = 
C\lo\ and therefore to an action of Eq. (|38[) which 
is non-local in time: 



S2 



drdr' 



ip(r) - <p(t') 



and is similar to the one of Caldeira and Leggett in the 
frame of quantum dissipation. This terminology might 
be misleading here, as the dissipation arising in the ef- 
fective action is not related to an external bath coupled 
to the system (phonons) . It is intrinsically related to the 
incommensurate CDW, i.e. the action is that of generic 
gapless modes. It may be understood by thinking about 
the CDW away from the impurity as an effective (inter- 
nal) bath governing the dynamics of the phase at the im- 
purity. Tracing out the internal degrees of freedom away 
from the impurity, while keeping the phase at the im- 
purity fixed, leads to non-trivial dynamics for the latter. 
This proce dure is well-known in ID disordered systems, 
cf. Refs. |5l l35l36| . the first extension to quasi-lD systems 
appearing in Ref. [2||. Because the phase is bounded 
along the time trajectory, i.e. Atp = 2ir, this part of the 
action is logarithmic in time: 



s 2 = b VYC Hl/uP), 



(40) 



where u = y/Y/C is the velocity of the collective electron 
structure and l/u an upper-time cut-off preventing the 
divergency of the action. 

With the help of Eqs. l|3*Tf and ijltty. minimizing the 
total action s — s\ + S2 with respect to (3, the optimal 
time reads: 



opt 



and the optimal action reads: 



r YC MlW b /bY). 



(41) 



(42) 



In Eq. I|42l) . we have assumed that the logarithm gives 
the major contribution with respect to unity. This is 
effectively the the inequality: Wbl/b ^> Y, is 

equivalent to the requirement that we are in the strong 
pinning regime. The latter also implies that the impu- 
rity strength and bulk modulus are determined by local 



electrostatics: Wb ~ e 2 / nb and Y ~ e 2 / nb, as may be 
check by dimensional arguments. The arguments of the 
logarithm therefore corresponds to: Wbl/bY = Ijb ^> 1. 

The second check deals with our initial assumption 
that at small times, the non-trivial dynamics of the phase 
at the impurity position are irrelevant. We may check 
that this is the case for an infinite, delta-function-like, 
back-scatterer (W& — > oo). In this case, the tunneling 
is instantaneous and the phase effectively jumps by 27T, 
i.e. by one period at the impurity position. The Ansatz: 
<p(r) = 2it9(t), where 6{t) is the Heaviside function, is 
then perfectly justified. This is also the case if the opti- 
mal time is smaller or equal to the smallest time-scale of 
the problem: /3 opt < a/u. Using again the fact that the 
plasmon velocity is related to our parameters with the 
help of: u = y/Y ~/C, this condition reads: Wb > Yb/a. 
The lower boundary is satisfied in the case of strong in- 
dividual pinning, W b ^ Y ^ e 2 /nb, for sufficiently dilute 
electronic systems (hence large r s , see below). 

Finally, we notice that the phase-phase correlation 
function in the absence of disorder scales as: < (p 2 >= 
h/b^YC (where h has been restored for clarity). The 
weakness of quantum fluctuations, on which our start- 
ing semi-classical approach was based, implies that: 
b\jYC h. This condition fulfills the requirement that, 
in Eq. (|42ll . the optimal action, s opt ^> 1. Because 
y/Y/C, and using again the fact that: Y ~ e 2 / ' Kb, 

2 /k/Hu = Uh/ef, where 



u 



we see also that: 



Uh = e 2 / na is the Coulomb energy scale and ep = Hu/a 
is the kinetic energy scale (a is the average distance be- 
tween electrons along a wire). We may therefore intro- 
duce the well-known parameter r s which is defined as: 



— =bVYC/h^> 1. 



(43) 



It follows from Eq. I|43|l . that the optimal action of 
Eq. ma Y be re-expressed as: 



"opt 



r s ln(Z/6), 



(44) 



where we have used the fact that Wb ~ Y and returned 
h to unity. It is crucial to notice that this action is non- 
WKB like because the strength of the barrier, Wb, is 
in the logarithm, as known from Larkin and Lee in the 
strictly ID case [34|]. The difference with the ID case 
here is that the inter-chain distance, b, appears in the 
logarithm instead of the average electron distance, a, in 
ID. From Eq. JSJl and |@3J|, this yields the localization 
length in the large impurity concentration case: 



1 



r s bN ln(l/M) 2 )' 



N-»N S 



(45) 



where N = I /lb has been used. 

We include now the long-range Coulomb. Following 
Eq. (|T2^ and the discussion below it, the inclusion of 
the Coulomb interaction amounts to replace: q 



~2 



q 2 x by 9 2 e 2 (q) 



q 2 + q 2 x /r D q -> g 2 [l + l/r D q], 
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q 2 . Because qrn <C 1 and at large 



where q 2 = q x , Hy 
distances along the chains, q x <C q y , this reduces to the 
non-analytic shift: q 2 — > q 2 /rD\q y \ in Eq. J3H|), which 
reads: 



dq 



1 



(2 7 r) 2 Co; 2 + r (7 2 /rD |^|- 



It is straightforward to show that the integration yields: 
J u = 'i-nb\/YC\uj\, which is the above result up to a nu- 
merical factor. In 2D and in the large impurity concen- 
tration case, the Coulomb interaction does not modify 
the result of Eq. (gSJ. 

The case where N <C N s requires stronger inter-wire 
couplings and will be considered next. 



8. The case of small impurity concentration 

In the small-impurity concentration case the inter- wire 
coupling is crucial, i.e. it gives rise to the non-trivial scale 
l s . For coupled wires, we are only aware of the results of 
Ref. j2|| dealing with the equivalent 3D geometry. In our 
2D case, including inter-chain interactions, and focusing 
on large-distances along the chains, q x <C \faq yi Eq. l|3l3|) 
reads: 



JZ 



rfq 



1 



This yields: 



Jrj — 



(2tt) 2 Cuj 2 



E s 



aYq 2 ' 



ln(E s /bVYC\w\) 



(46) 



where E s oc e 2 / kI s is the soliton energy. The total action, 
s = si + S2, in the logarithmic approximation, therefore 
reads: 



W b (3- 



E s (3 



\n(E s (3/bVYC) 



The logarithmic factor provides a minimum to the action 
which, in the strong pinning regime, Wb S> E s , is then 
given by: 



^opt 



Wb 
E/ 



(47) 



where Eq. i|43|) has been used. 

Next, we include the effect of the Coulomb interaction 



shifting q 2 to g 2 62(q) 



q 2 



ql/roq, cf. Eq. {T^J and 



the discussion following it. At large distances along the 
chains, q x <C \/aq y , the shift reduces to q 2 — > aq 2 + 
q 2 /r D \q v \. The kernel of Eq. then reads: 



JZ 



(iq 



(2tt) 2 Cuj 2 
The integrations yield: 



aYq 2 + Yq 2 /r D \q y \ 



E» 



1 - {by/YC\u\/E a y/* 



(48) 



Eq. (|48(l shows that, in the presence of the Coulomb in- 
teraction, the dynamics enter the kernel in a perturba- 
tive way with respect to the static part E s . This was 
not the case in ID as well as in 2D without the long- 
range Coulomb interaction, cf. the logarithmic factor in 
Eq. H46JI . The expansion of the kernel, up to second or- 
der, leads to the following total action: 



!YC \n(l,/u0), 



where our cut-off is now l s . The optimal time is found to 
be: 



opt 



W b 



(49) 



which leads to the following optimal action: 



5 opt 



In the strong pinning regime, Wb ~ Y, and with l s ^> b, 
the action is dominated by the logarithmic contribution 
which, up to a numerical factor, reads: 




^opt 



r s \n(l s /b). 



(50) 



where Eq. (|4^)l has been used as well as the strong pin- 
ning result: Wb ~ Y . The long-range Coulomb therefore 
returns us to the logarithmic action which is non WKB- 
like. 

From Eqs. I|35() and l|50[) . the localization length in the 
small impurity concentration case reads: 



i 



N < N, 



(51) 



r s WV s ln(l/iV s 6 2 ) : 
where N = 1/lb has been used (N s — l/l s b). 

9. General expression and remark 

Eqs. Q45f) and (|51|l yield the general expression for the 
localization length: 



£~ x = r s b max{A^,iV s } ln(- 



1 



t{N,N s }b 2 



(52) 



which we will use in the following section. 

As a final note, we follow Ref. 26], where it has been 
mentioned that the collective dynamics of the CDW may 
include a contribution from the amplitude, |A(r)|. This 
contribution arises because the order parameter is com- 
plex: A(r,i) = \ A\exp(i(p). At zero temperature the am- 
plitude mode is frozen but in VRH we may be interested 
in reaching the thermally activated regime. In this case, 
there is an additional regular kinetic energy contribution 
oc Iuj 2 , where the momentum of inertia / depends on fluc- 
tuations of the amplitude. This brings an additional con- 
tribution to the total action of Eq. (j5U|) : oc 1/(3. Eq. lf5U|) 
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would then be valid provided that: / <C b 2 YC jWb- In where the parameter Tm reads: 



the other case, I ^> b 2 YC /Wt, the final action reads: 



'opt 



2\/lW~ b 



(53) 



where the second term is a correction and L — min{Z, l s }. 
Notice that the first term, with a square-root dependence 
on the barrier, is WKB-like. In the following, we will 
assume that the amplitude mode has a negligeable con- 
tribution and we will use the final expression for the lo- 
calization length, Eq. I|52|l . based on the non WKB-like 
actions. Depending on the system under consideration 
the WKB contribution may however play a significant 
role. One would then have to use Eq. (|53|l together with 
Eqs. (|34fl and l|35|> to derive the localization length. 



V. LINEAR VRH LAWS 

This section deals with the hopping laws of the quasi- 
1D Mott- Anderson insulators in the linear response 
regime. Our results depend crucially on the Coulomb 
gap shapes derived in Sec. IIIII and on the impurity- 
dependence of the longitudinal localization length de- 
rived in Sec. IIVI Notice that all transport laws below, 
and especially the parameter characteristic of these laws, 
are given up to a numerical coefficient of the order of 
unity. 

Recall that the semi-phenomenological arguments ini- 
tially introduced by Mott in order to derive the dc con- 
ductivity, <r(T), are based on minimizing (with respect 
to the coordinates x and y) the toy-action: 



S = 2x/£ x + 2y/£ y + T/T, 



(54) 



where the quantum and classical parts are related by 
a Fermi's Golden rule involving the density of localized 
states: 



vTS w 1, 



(55) 



where S — xy is the area. Eqs. (|54() and l|55[) describe 
the process with the help of which an electron may hop 
from a site i to the site j within a disorder band with 
a constant density of localized states, v. The first two 
terms in Eq. I|54|) represent the overlap between the two 
states separated by x and y. Because the system is dis- 
ordered the energies of these states are different and a 
phonon has to be involved in the process of hopping. 
The thermal energy, T, is then determined with the help 
of Eq. I|55|l and corresponds to the average energy spac- 
ing between the localized states involved in the hopping: 
r = 1/vS = 1/vxy. Introducing this value of T in 
Eq. (|54() and minimizing the total action with respect 
to x and y yields the optimal hopping distances: 

x opt = C (TAf/T) 1 / 3 , y opt = £ y (Tm/T) 1 / 3 , 



(56) 



This yields an optimal action: S op t oc {Tm/T) 1 / 3 . The 
corresponding d.c. conductivity a cx exp(— S op t), is the 
Mott law for variable-range hopping and reads: 



<tm(T) = a (r)exp[-(T M /T) 1 / 3 ] 



(57) 



where ctq is a temperature-dependent pre-factor (the 
temperature-dependence arising from the phonon- 
scattering time which is a power- law). This law may 
be straightforwardly extended to any dimension. It may 
have a disorder-dependence through the constant DOS 
of the disorder band, v. 



N K 
e 2 /nl e 2 b' 
N k N 



V cx 



e 2 /nl s e 2 b N„' 



N^>N S , (58a) 
iV<A s . (58b) 



This estimate was derived on the basis that we have N 
localized states per unit volume, each with an energy 
e 2 / 1 kI in the limit N 3> N s (charge excitations extend 
over segments between impurities, I) and e 2 / kI s in the 
limit N <C N s (charge excitations have their length-scale, 
l s ). This shows that v grows with N in the small impurity 
concentration case and saturates when N becomes larger 
than 7V S . An additional dependence may come from £ x . 
In particular from Eqs. H5(jfl and l|52l) . the parameter of 
the Mott-law for 2D quasi-lD systems reads: 

r s bNln(l/Nb 2 ) , s 

T M = H ~, N^N S , (59a) 



T, 



M 



r s bN s ln(l/7V s 6 2 ) 



, N < N s , (59b) 



where N = 1/lb is the impurity-density, b the inter- wire 
distance and JV S has been defined by Eqs. I|1U|) and i|17fl . 
It is interesting to notice, from Eq. (|59a|l . that in the large 
impurity concentration case there is a linear dependence 
on N arising from £ x [y is constant from Eq. I|58a|) ). 
On the other hand, in the small impurity concentration 
case, Eq. (|59b() . £ x saturates and we have Tm ~ l/N 
from the N— dependence of v, Eq. (|58b(l . In this case, 
the conductivity increases with increasing disorder. 

When the long-range Coulomb interaction is taken into 
account, the latter leads to a depletion of low-energy 
states in v. As we have seen in Sec. IIIII this depletion 
corresponds to a soft Coulomb gap g(e) in v. Following 
Efros and Shklovskii, in the presence of this Coulomb 
gap, one has to replace Eq. by: 



ff(T)rs«i. 



(60) 



With the help of Eq. I|25|l , corresponding to a usual ES 
Coulomb gap shape, and the arguments above, the d.c. 
conductivity is the Efros-Shklovskii law for VRH: 



&es{T) = a (T)exp[-(T £S /T 1 ) 1/2 ] 



(61) 



12 



where the parameter Tes depends on the specific 
Coulomb interaction among charge carriers: 



Tes — 



K \/ £x£,y 



(62) 



Substituting the expression of the localization length, 
Eq. this parameter reads: 



T, 



ES = — 1 

K 



lr s Nbln(l/Nb 2 ) 



Tes — 



lr s N s b\n(l/N s b 2 ) 



N > N s , (63a) 



N < N s . (63b) 



In the large impurity concentration case, N 3> N s , this 
parameter increases with the impurity concentration, N, 
which leads to a decreasing conductivity as a function of 
disorder. This impurity-dependence originates from £ x . 
In the small impurity concentration case, N <C N s , £ x 
saturates and TeSi as we U as the corresponding conduc- 
tivity, become disorder-independent. 

At higher energies, the usual Coulomb gap of Eq. i|25|) 
crosses over to the unusual Coulomb gap of Eq. (|2*7Jl . The 
latter yields a conductivity that we denote as: oi, and 
reads: 



C r 1 (T) = C r (T)ex P [-(T 1 /T) 3 / 5 ] 



(64) 



The unusual exponent, i.e. 3/5 instead of 1/2, is related 
to the unusual Coulomb gap of Eq. 127|) . Furthermore, 
the parameter T\ reads: 



Ti 



(bn 2 K x £,x£,y) ^ 



(65) 



Substituting the expression of £ x , Eq. l(52*|l . and of k x , 
Eq. yields: 



Ti = — - Nb 2 
Kb 



Ti = — 

K 



r s b\n{l/Nb 2 ) 
Tv 



r s bN 2 N s ln{l/N s b 2 ) 



1/3 



iV> iV s ,(66a) 



1/3 



, N < 7V s .(66b) 



The Ti parameter depends on disorder through both £ x 
and Ka,, in the large- impurity concentration regime. On 
the other hand, £ x saturates in the small impurity con- 
centration case and the disorder dependence of T\ orig- 
inates then only from k x . In both cases this disorder- 
dependence of Ti yields a conductivity, <ji, which de- 
creases with increasing disorder. 

We then focus on the peculiar 2/3 Coulomb-gap sector 
in the small impurity case, cf. Eq. I|31|) . This Coulomb- 
gap yields: 

a a (T) = a (T)exp[-(T 2 /T) 6 ' u ], N « N s . (67) 

The unusual exponent, i.e. 5/11 instead of 1/2, is again 
related to the unusual Coulomb gap (here oc e 2 / 3 ) of 



T 

T, 

i) 



Activated 




1/JV, 



1/N 



FIG. 4: The generic phase diagram of a 2D quasi-lD system 
as a function of temperature T and the inverse average im- 
purity concentration 1/N [T s oc E s the soliton energy]. This 
phase diagram displays the VRH laws with exponents 1/2, 
3/5 and 5/11 as well as Mott law and activation at higher 
temperatures. 



Eq. I|31|) . Furthermore, the parameter T 2 reads: 

"I 1/5 



To 



ar D 



£3 £3 



N<^N S 



(68) 



Substituting the expression for £ x> Eq. (|52|l . this param- 
eter reads: 



To 



ar 3 s b 4 N 3 In 3 (1/N s b 2 ) 



,1/5 



(69) 



This parameter, as well as the corresponding conductiv- 
ity, (72, are disorder- independent. 

Finally, at high energies (within the constant DOS 
v), these VRH laws cross-over to the so called nearest- 
neighbor hopping (NNH) or activated law, see Ref. for 
a review. The latter reads: 



(?nnh(T) = o- (T)exp[-E NNH /T], 
where the NNH energy reads: 



Ennh 
Ennh 



e 2 bN/n, N > N s , 
e 2 bN s /n, N <C N s . 



(70) 



(71a) 
(71b) 



These results are summarized on the temperature — 
impurity-concentration phase-diagram of Fig. The 
crossover lines between the different laws in Fig. 0] are 
determined with the help of the following arguments. 

The expression of ei, Eq. J2HJ, determines the 
crossover temperature lines Ti(iV) between the ES and 
the 3/5-laws [e x = T^Tes/T^) 1 ' 2 ] with Tes given by 
Eq. H63|) ■ This crossover line reads 



T cl 



nb 



(Nb 



2.7/2 



t c1 = - 2 (my 



,bN s log( 



l 



iV> iV s ,(72a) 
-, 7V< iV s (72b) 



In the case of large impurity concentration, the expres- 
sion of Ai, Eq. H28fl determines the total width of the 
Coulomb gap and depends on the constant DOS v. This 
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width determines the high temperature crossover line 
T c2 {N), [Ai = between the 3/5-law and 

nearest- neighbor- hopping, cf. Eq. (|71a|l . for N ^>> N s . 
With the help of Eq. JH1, T c2 (N) reads: 



T c2 (N)= -Nb* 
A> A, 



5/4 



r s 61n(l/A& 2 ) 



(73) 



Using the estimation for v, Eq. I|58af) . Eq. l|73Jl may be 

re-expressed as: 



T c2 (N) = — Nb 2 

KO 



r s bln(l/Nb 2 ) 



N->N S . (74) 



From Eqs. (|73a|) and (73J, we see that: T c2 (N)/T cl (N) = 
1/(6 2 A) 5 / 2 > 1, so that at the boundary N ~ N s , the 
corresponding crossover temperatures are separated by 
the large dimensionless factor (Z s /6) 5 / 2 = 1/a 5 / 4 . More- 
over, we define T s = T c2 (N s ) which reads, up to a loga- 
rithmic factor: 



T s = T c2 (N s ) (xE s , 



(75) 



where: E s = e 2 /nl Sl is the creation energy of a soliton. 
With these notations we have: 



T cl (A s )cxa 5 / 4 T s «T s 



(76) 



We could also show that in this large impurity concentra- 
tion case there is no room for the Mott law of Eq. (|59a|l . 
i.e. the crossover line from the 3/5— law to the Mott law 
coincides with the crossover line from the Mott law to 
the activated- or NNH- law. 

In the small impurity concentration case we have al- 
ready given the expression for the T\— line in Eq. (|72b() . 
which matches smoothly Eq. I|72a|) at N ~ N s . On 
the other hand, the T 2 — line splits into three crossover 
lines thereby opening two sectors for the new 5/11— law 
and the Mott law. The energy e 2 , Eq. I]32p. deter- 
mines the temperature crossover-line: T® 2 (N), [e 2 = 
T c2 (Ti/r c 2) 3/5 ], between the 3/5-law and the 5/11-laws 
for N -C N s . With the help of Eq. l|66rJ)> for the param- 
eter Ti and Eq. T? 2 (N) reads: 



N 



A s 



11/4 



r s 61og(l/A s & 2 )' 



A< A s 



(77) 

which matches smoothly Eq. (|74|l at N ~ N s and de- 
creases more abruptly with 1 /N for N -C N s . 

At the next crossover-line: 7j 2 (iV), the 5/11— law 
crosses over to the Mott VRH law with the parameter 
Eq. l|59bl) . This crossover-line reads: 



T; 



hi o 



N 



5/2 



r s blog{l/N s b 2 )- 



A< A s , 
(78) 



which matches smoothly Eqs. (|74|l and l|77|) at N ~ N s . 
Finally, at higher temperatures a crossover-line bridges 
the Mott law of Eq. (pp with the NNH law of Eq. ffflfo : 



T C C 2 (N) 



Kb a 



N 

n1 



1/2 



r s 61og(l/A s 6 2 )' 



A< A s 



(79) 



which matches smoothly Eqs. ifTIj). f77|l and (JJSJ at AT . 



VI. NON-LINEAR VRH LAWS 

This section deals with the hopping laws of the quasi- 
1D Mott-Anderson insulators in the non-linear response 
regime. Our results depend crucially on the Coulomb 
gap shapes derived in Sec. IIIII and on the impurity- 
dependence of the longitudinal localization length de- 
rived in Sec. IIVI Our arguments follow closely those 
of the linear-response regime of Sec. El Notice that all 
transport laws below, and especially the parameter char- 
acteristic of these laws, are given up to a numerical coef- 
ficient of the order of unity. Notice also that our present 
arguments do not allow us to determine the pre-factor of 
the current; we therefore focus only on its main exponen- 
tial dependence. 

The VRH laws derived in Sec. are valid in the linear 
response (or ohmic) regime where the current is linear in 
the applied electric field, j — a(T)£. At a given tempera- 
ture, upon increasing the electric field, a crossover should 
take place, below the threshold for the sliding of the elec- 
tronic crystal, to a non-linear regime. In the frame of 
doped semiconductors, such a transition has been stud- 
ied in Ref. |3?] by extending the Mott argument presented 
in the previous subsection on the linear VRH laws. Such 
an extension amounts to replace, in Fermi's Golden rule 
of Eq. H55f) . the hopping energy T by the energy provided 
by the electric field, e£r, during the motion of an elec- 
tron along a distance r from the initial localized state. In 
d dimensions, the optimal hopping distance is therefore 
given by: r opt = \/{ve£) 1 l 1+d . Substituting this value 
in the tunneling probability oc exp(— 2r/£) yields: 



3ms(£) ~ exp[ 



/ £ Al S \ — - — i 
-( ) 1 + d 



£ 



1 



MS 



l+a 



(80) 



where the index MS refers to Mott-Shklovskii, see 
also Ref. IsH- The crossover between the linear and 
non-linear regimes takes place when: (£ms/ £) 1 ' 1+d = 
(T M /T) 1 / 1 + rf , where the parameter of the Mott law: 
Tm = l/^C d > m d-dimcnsions. At a given temperature, 
the threshold field for the non-ohmic regime is given by: 
£ c = (£m/Tm)T, that is: 



e ^ 



(81) 



in all dimensions. This returns us to the arguments of 
the Introduction, see Eq. J3J and discussion around it. 
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It is straightforward to generalize such results in the 
presence of a Coulomb gap. We need simply to replace 
the constant DOS, u, in Fermi's Golden rule by the cor- 
responding Coulomb gap: g(s), evaluated at the energy 
of the charge excitation: e = eEr. The generic form of 
the Coulomb-gap (CG) in isotropic system is given by 
Eq. 124(1 . This leads to an optimal hopping distance in- 
dependent of dimensionality: r opt = (e 2 / KeE) 1 ! 2 . The 
corresponding current is therefore given by: 



j (8) ~ exp[-(£ /£) 2 



So 



net; 2 



(82) 



Strictly speaking, Eq. (|S2)l is valid for d > 1. In 
d = 1, there is a logarithmic Coulomb gap: g(e) — 
v I log(e 2 iV6/re|e|), see Ref. Q for a similar derivation 
of this law in the ohmic regime. For the ID case, the 
parameter of Eq. I|82|) is therefore given by: 



So 



log[(e 2 M>/ K ), 



vet, 1 



(83) 



in the logarithmic approximation where the energy of the 
charge excitation has been taken equal to e — e£r opt and 



r opt = \/\/ve£ is the optimal hopping distance in ID. 

These results may be extended to quasi- ID systems. 
For simplicity, we assume that the electric field is par- 
allel to the chains. More general results with the two 
components of the field can be derived in the same way. 
We also focus on the 2D case. Then, Fermi's Golden 
rule of Eq. 1(55(1 is generalized to: v(e£ x x)xy ~ 1. The 
tunneling action is anisotropic and takes the usual form 
S = 2x/£ x + 2y/£ y . Substituting the a:— component of 
the hopping length in this action and minimizing the re- 
sulting expression with respect to y yields to the optimal 
hopping distances. Substituting the latter in the action 
we finally obtain a current: 



3ms{S) ~ exp[-(— -)*], £ MS = c2c ■ 



(84) 



Substituting the expression for the longitudinal localiza- 
tion length, Eq. i|52|). this yields: 

£ms ^ fc W(i/iv^ w>>Wji 



Sms = 



r 2 s b 2 N 2 \n 2 {l/N s b 2 ) 



N <iV s 



Substituting the expressions for u, Eqs. (JSSJl, yields the 
explicit N— dependence of the parameter: 



-MS 



Kb 



(iv& 2 ) 2 -2 — 



N^N S , (86a) 



e 2 N< r 2 (N s b 2 ) 2 In 2 
Kb iv et, y 

The next step consists, as for the isotropic system, in 
introducing the Coulomb interaction, which amounts to 



replace the constant DOS, v, by the Coulomb gaps deter- 
mined in Sec. 11111 For the linear Coulomb gap of Eq. i(25|) 
the non-linear law follows Eq. ((82(1 with a modified pa- 
rameter: 



1 



(87) 



Substituting the expression for the longitudinal localiza- 
tion length, Eq. (|52() . this parameter reads: 



£o = 



£o= P - 



r s bN ]xi(l/Nb 2 ) 
7173 

52/ 

r s bN s ln(l/iV s 6 2 ) 
?3/ 



3/2 



, N > N s , (88a) 
3/2 

, N <c iV s . (88b) 



On the other hand, the anomalous quadratic CG of 
Eq. ((27(1 leads to a 3/5— law in the non- linear regime: 



h{S)^eM-{S 1 /Sf'% 
where the parameter reads: 

1 1/3 

K 



K x£, x £yb 



(89) 



(90) 



Substituting the expression for the longitudinal localiza- 
tion length, Eq. (|52() . this parameter reads: 



2 ,2 



£ = ^ (iViV 2 ^ 



N^N S , (91a) 



,7V<^ s .(91b) 



Finally, in the small impurity concentration case, the 
2/3— Coulomb gap of Eq. 1)31(1 leads to a 5/11— law for 
the current: 

j 2 (£) ~ exp[-(£ 2 /£) 5 / u ], N^N S , (92) 
with a parameter 

"I 1/5 



£2 = 



ar D 



£8 £3 



, N<£N S 



(93) 



Substituting the expression for the longitudinal localiza- 
tion length, Eq. ((52(1 . in the small impurity case, this 
parameter reads: 



£2 



kI 2 



fr 8 s b 3 l n 8 (l /N s b 2 ) 
^y 



1/5 



N<^N S 



(94) 



From Eqs. (fHH|> . (JSSJ, (jHU and (|M|l for the parameters 
and the expressions of the corresponding currents, the 
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In a 




5/11 Mott Act. 

1/A r 



FIG. 5: The logarithm of the resistivity, in the linear re- 
sponse regime, versus the inverse average impurity concen- 
tration, 1/N, for a 2D quasi-lD system at a temperature: 
T < q 5/4 T s , where T s <x E s , see Fig. The ES [Eq. (STJ], 
3/5 [Eq. fig}], 5/11 [Eq. {§3], Mott [Eq. igZJ] and activation 
[Eq. 17UH 1 law succeed each-other with decreasing N. 



non-monotonous behavior of j as a function of N is rather 
clear. 

We shall not detail the field crossover lines between the 
various laws. This can be done exactly in the same way as 
for the linear-response regime case. Rather than that we 
would like to point out that the phase-diagram, electric- 
field vs. impurity concentration, may be obtained from 
the one of the linear law by assuming that the electric 
field gives rise to an "effective" temperature, Tg. This 
effective temperature reads: 



T E = eS 



'opt ) 



(95) 



where x op t is the optimal hopping distance along a wire 
(recall that the electric field is assumed to be parallel to 
the chains). For a given non- linear (NL) law of parameter 
£nl and exponent Jnl, the optimal hopping distance 
along the wires is defined as: 



(N,S) = ( X (£ NL /£) 



1NL 



(96) 



and depends on the applied electric field and, eventu- 
ally, the impurity concentration. The phase diagram, 
Te(N) vs. N, therefore has a non-trivial re-scaling via 
the TV— dependence of Tg. 



VII. CONCLUSION 

We attempted to construct a semi-phenomenological 
theory of variable-range hopping for 2D quasi- ID systems 
such as arrays of quantum wires in the Wigner-crystal 
regime. We have closely followed the phenomenological 
arguments of Mott, Efros and Shklovskii to derive the 
Coulomb gap shapes, Sec. IIIII as well as the main expo- 
nential dependence of the transport laws in the linear, 
Sec. [V] and non-linear, Sec. I VII response regimes. Our 
approach has been supplemented with some microscopic 
arguments necessary to derive the impurity-dependence 
of the longitudinal localization length, Sec. IIVI Both 
Coulomb gap shapes and transport laws were found to 
have rather unusual features with respect to known re- 
sults in the field of disordered semiconductors 91. These 



unusual features arise because of the non-trivial dielec- 
tric properties of the systems under consideration, SeclTTl 
They are two-fold: a non-monotonous dependence of the 
conductivity or current as a function of disorder and a 
highly non-universal exponent 7. In the linear response 
regime, the richness of exponents is displayed on the 
phase diagram of Fig.QJ Despite the fact that some expo- 
nents may be close to each-other (1/2 = 0.50, 3/5 = 0.60, 
5/11 = 0.45, 1/3 = 0.33), the corresponding monotonic- 
ity of the conductivity as a function of disorder enables 
further discrimination between the various laws. This is 
schematically displayed on Fig. El Moreover, 2D layers of 
wire arrays are experimentally accessible, to our knowl- 
edge 15], but we are unaware of systematic transport 
experiments in the strongly localized regime for such sys- 
tems. Our theory is therefore only of predictive nature. 
We hope, however, that it will be of some interest to both 
theorists and experimentalists working in the field. 



APPENDIX A: 2D MODEL OF INTERRUPTED 
STRANDS 

As an alternative to the use of Eqs. l(T%|) and ijT§|) 

the "disorder averaged" electrostatic potential may be 
derived following the results of the interrupted strand 
model, see Ref. 1311 of Rice and Bernasconi for the 3D case. 
In this approximation, we assume that the dielectric con- 
stant along the chains is given by Eq. 11221) . This equation 
shows that for L x /ro = l/ro 3> 1, the contribution of 
the electronic part along the chains is much larger than 
the transverse part as well as the host lattice dielectric 
constant (n y — k <ti k x ). The problem therefore reduces 
to determine the electrostatic potential of a charge car- 
rier in a layer of a quasi-one dimensional system whose 
longitudinal dielectric constant is much larger than the 
dielectric constant of the surrounding media. Keldysh 39 
has solved a similar problem for an isotropic layer. In 
the present case, the Fourier transform of the interaction 
potential at (r, z) due to a point charge at the origin of 
the layer is found by usual methods of electrostatics and 
reads: 



Aire 2 cosh (#6/2 + S) 
q + q smh(qb + 25) 



exp(q(b/2-z) + S), 
(Al) 



where q = (q x , q y ) is the two-dimensional reciprocal vec- 
tor of r, q 2 = n x q 2 + q 2 and 5 is given by: 



log 



(A2) 



At large distances qb <C 1, i.e. for x ^> I, that will be of 
interest to us in the following, and for z = 0, Eq. (|A1(I 
reduces to Eq. l(T%|l with the following dielectric function: 



e 2 (q) « k[1 



Kx bq 2 x . 
K q 



(A3) 
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This expression is equal to Eq. iJTSj l in the limit q = q x <C 
L~\ with k x ss K,{L x /r£)) 2 and rry w b. This shows that 
in this limit, both approaches are equivalent. 
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